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Abstract 

We consider the interplay of fundamental and matter-induced T violation effects 
in neutrino oscillations in matter. After discussing the general features of these effects 
we derive a simple approximate analytic expression for the T-violating probability 
asymmetry AP^ for three-flavour neutrino oscillations in a matter with an arbitrary 
density profile in terms of the two- flavour neutrino amplitudes. Explicit examples are 
given for the cases of a two-layer medium and for the adiabatic limit in the general case. 
We then discuss implications of the obtained results for long baseline experiments. We 
show, in particular, that asymmetric matter effects cannot hinder the determination 
of the fundamental CP and T-violating phase 5c p in the long baseline experiments 
as far as the error in this determination is larger than 1% at 99% C.L. Since there 
are no T-violating effects in the two-flavour case, and in the limits of vanishing #13 
or Amj] the three-flavour neutrino oscillations effectively reduce to the two-flavour 
ones, studying the T-violating asymmetries AP^ can in principle provide us with a 
complementary means of measuring #13 and Am^. 
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1 Introduction 



The experimental evidence for neutrino oscillations implies that neutrinos are massive and 
that mixing angles and CP phases exist in the lepton sector. The measurement of these 
neutrino parameters is very important, not only because it will provide us with an informa- 
tion about neutrino properties, but also because it may have interesting implications for the 
structure of theories at very high energies, which might "explain" the low energy parame- 
ters. Since neutrinos could in general have both Dirac and Majorana masses, they might, 
in principle, also allow more insight into the flavour problem than quarks do. Extraction 
strategies for neutrino masses and mixings in current and future experiments are therefore 
an important subject under study. Future experiments also offer to study CP, T, and CPT 
properties in the neutrino sector. Local quantum field theories are, however, for general 
reasons invariant under CPT, and any CP violation implies then a correlated T violation. 
An independent study of T violation might therefore not appear very interesting, unless 
it is sensitive to tiny CPT violating effects induced at very high energy scales by physics 
going beyond local quantum field theory. In this paper, we discuss that this is different for 
neutrinos and that an independent study of T violation for neutrino oscillations in matter 
offers interesting insights. The point is that the presence of matter in the experimental 
setup violates by itself CP as well as CPT and gives thus rise to an extra CP violation in 
addition to the intrinsic CP violation. Furthermore, there are extra T-violating effects if 
the matter density profile seen by neutrinos is asymmetric. This extra T violation does not 
follow directly from CPT because CPT itself is violated by the presence of matter. 

One can understand these effects by comparing CP, T, and CPT properties of neutrino 
oscillations in vacuum and in matter. Using CPT one finds in vacuum P[y a — > Vb) = 
P{vb — > v a ). This is no longer true when one studies neutrino oscillations in matter, i.e., a 
CP- and CPT- asymmetric environment, where the oscillation probabilities of neutrinos and 
antineutrinos change in a different way due to the differences in coherent forward scattering 
in a given medium (MSW effect). As a consequence, CPT is violated by matter effects, 
i.e., P{y a — > Vb) 7^ P{pb v a \ Moreover, the total CP violation is now a combined effect, 
where intrinsic CP violation must be separated in analyses from effects induced by the 
CP- violating environment. A quantity which measures the total CP violation is AP^ b p = 
P{y a Vb) — P{p a —> Vb). Such "CP differences" or corresponding asymmetries have been 
studied for appearance channels like v e — > at future long baseline experiments (see, e.g., 
Refs. |j]-f|]), and one can nicely see how the CP difference is in general a combination of 
intrinsic CP effects and matter effects, while in vacuum AP^ P depends only on intrinsic 
CP violation. Note that one can easily see from CPT and the definition of AP^ P that in 
vacuum CP violation can only occur in appearance channels, i.e., for a ^ b, while in matter 
one has in general AP^ a p ^ 0. 

In analogy to the CP difference AP^ b p one can also define a "T difference" APj b = 
P(v a -> Vb) ~ P(v b -> v a ) and the "CPT difference" APg PT = P(u a -> v h ) - P(u b -> v a ). 
In vacuum, where CPT holds, one has AP^ P + AP^ = AP^ b PT = and AP^ is thus given 
directly by AP^ b p . Since in matter CPT is no longer valid, one now has AP^ b PT ^ 0. For 
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this reason, T violation is no longer correlated with CP violation. In particular, matter does 
not necessarily break T invariance; as was mentioned above, this only happens if its density 
profile seen by neutrinos is asymmetric. 

T violation in neutrino oscillations in matter has been studied in a number of papers 
(see, e.g., Refs. P-pO|). However, in most of these studies matter of constant density was 
considered, where only intrinsic T violation is possible. T violation in neutrino oscillations in 
non-uniform matter was considered in Refs. [^,^0[. In the first of these papers T- violating 
effects of solar matter on the oscillations of O(GeV) neutrinos, which could be produced 
in the annihilation of WIMPs inside the sun, are considered, whereas in the second paper 
small perturbations of constant matter density profiles were discussed. In the present paper 
we study T violation in neutrino oscillations in matter of an arbitrary density profile and 
discuss the interplay of the fundamental and matter-induced T violation. In particular, we 
discuss where the matter-induced T violation effects may play a role and where they can be 
safely ignored. 

The paper is organized as follows. In Sec. ^| we discuss the general features of T viola- 
tion in neutrino oscillations in matter. Subsequently, we give approximate analytic results 
for arbitrary matter density profiles in Sec. |3] and discuss the implications for long base- 
line experiments, solar, atmospheric, supernova, and cosmological neutrinos in Sec. f|. We 
discuss the obtained results and conclude in Sec. |. In Appendix A we give details of our 
general analytic approach, whereas in Appendix B its application to two particular cases is 
considered. 

2 General features of T violation in neutrino oscilla- 
tions in matter 

Oscillations of neutrino flavour in vacuum or in matter are described by the Schrodinger 
equation 



where \u) is the neutrino vector of state and H(t) is the effective Hamiltonian which in 
general depends on time t through the t-dependence [] of the matter density N(t): H(t) = 



We will assume that neutrinos are stable and are not absorbed in matter; in this case 
the Hamiltonian H(t) is Hermitian. It is convenient to define the evolution matrix S(t,to): 



: We consider the evolution of the neutrino system in time rather than in space (since for relativistic 
neutrinos r — t, both descriptions are equivalent). The coordinate dependence of the matter density profile 
N(r) therefore translates into its time dependence N(t). 

2 Neutrino flavour evolution in matter depends on the total matter density as well as on the electron 
number fraction Y e . Throughout this paper we will for brevity call them collectively "matter density 
profile" . 



i-\v) = H{t)\v) , 



(1) 



H[N(t)} 5 




(2) 



2 



It has the obvious properties 

S(t, t ) = S(t, tJSfa, t ) , S(t , t ) = l, S(t, t )S(t, t o y = 1 , (3) 

where the last property (unitarity) follows from the hermiticity of H(t). From Eqs. ([D and 
(H) it follows that S(t, to) satisfies the equation 

ij t S(t,t )=H{t)S(t,t ). (4) 

It is also sometimes useful to consider the evolution equation of S(t, to) with respect to its 
second argument, t . Differentiating the equality S(t,to)S(to,t) = 1 with respect to t and 
using Eq. (^), one finds 

i^-S(t,t ) = -S(t,t )H(t ) . (5) 

CLto 

The amplitudes of neutrino flavour transitions are just the elements of the evolution matrix 
S: 

A[v a (t ) -> Vb(t)] = [S(t,t )] ba - (6) 

We are interested in the properties of the solutions of the evolution equations ([!]) or (^) 
with respect to the time reversal transformation T. Unlike in the case of CP transformation, 
these properties cannot be directly experimentally tested as one cannot change the direction 
of time. However, time reversal can be studied by simply interchanging the initial and final 
neutrino flavours. This can be readily seen in the case of neutrino oscillations in vacuum, 
for which 

[S{t,to)] ba = Y J Ua ] e-^ t ~^W bj} (7) 

3 

where U and Ei are the lepton mixing matrix and energy eigenvalues in vacuum, respectively. 
Indeed, time reversal interchanges the initial and final times, to and t, respectively, which 
leads to the complex conjugation of the exponential factors in Eq. (|7|). On the other hand, 
the interchange of the initial and final neutrino flavours a <-» b (i.e., S —>■ S T , where the 
superscript T denotes transposition) is equivalent to the complex conjugation of the matrix 
elements of the leptonic mixing matrix U in Eq. (|7]). The action of these two operations on 
the elements of the evolution matrix S in Eq. (0) differ from each other only by complex 
conjugation. Since the transition probabilities P a b = P{y a — > V\,) are just the squares of 
the moduli of the matrix elements Sb a , these two procedures are physically equivalent. This 
means, as is well known, that in vacuum, instead of studying neutrino oscillations "backward 
in time" , one can study the oscillations forward in time, but with the initial and final flavours 
interchanged. 

The situation is, however, less obvious in the case of neutrino oscillations in matter 
of non-constant (and, in general, asymmetric) matter density. Time reversal means, in 
particular, that the matter density profile has to be traversed in the opposite direction. That 
is to say, one has to consider the evolution in the "reverse" profile N(t) rather than in the 
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"direct" profile N(t). Here N(t) should be understood as the profile seen by neutrinos when 
the positions of the neutrino source and detector are interchanged. Therefore, the following 
question arises: When one reduces the probability of the time-reversed neutrino oscillations 
in matter to that for the interchanged neutrino flavours, does one have to simultaneously 
replace the direct profile by the reverse one? There is some confusion in the literature 
regarding this issue, and therefore we believe that it is worth clarifying it here. In general, in 
the case of neutrino oscillations in matter of non-constant density, no closed-form expressions 
for the transition probabilities exist, and the arguments similar to those applied to Eq. ([?]) 
cannot be used. There is, however, a very simple and general argument which does not 
depend on whether neutrinos oscillate in vacuum or in matter and what the matter density 
profile is. 

Indeed, under time reversal the arguments of the evolution matrix S(t, to) are inter- 
changed: 

T: S{t,t ) - S(t ,t) . 

From Eq. (]3|) one then finds 

S{t , t) = S(t, to)' 1 = S(t, to) 1 " = [S(t, t ) T Y . (8) 

Therefore, the operations of time reversal and interchange of initial and final neutrino 
flavours [i.e., transposition of S(t, t )} are related by complex conjugation and so lead to the 
same transition probabilities in matter with an arbitrary density profile as well as in vac- 
uum. This, in particular, means that instead of considering neutrino oscillations in matter 
"backward in time" , one can consider the oscillations between interchanged initial and final 
neutrino flavours forward in time and in the same (i.e., direct) matter density profile. 

In vacuum, due to CPT invariance, T violation is equivalent to CP violation. In partic- 
ular, the CP-odd and T-odd differences of neutrino oscillation probabilities 

Apc p = P ab - P nl and APj 6 = P ab - P ba 

are equal to each other. In the case of n > 3 neutrino species they are expressed through 
(n — l)(n— 2)/2 Dirac-type CP- violating phases {Sep} 0- 

The situation with neutrino oscillations in matter is drastically different. Ordinary 
matter is both CP and CPT-asymmetric as it consists of particles (electrons and nucleons) 
and not of their antiparticles or, in general, of unequal numbers of particles and antiparticles. 
The violation of CPT by matter implies that CP and T violation effects are in general 
different. CP is violated by the very existence of matter. This violation manifests itself 
even in the two-flavour neutrino systems (or, in general, even in the absence of fundamental 
CP violation, i.e., when 5cp = 0) - matter may enhance the oscillations between neutrinos 
and suppress those between antineutrinos, or vice versa jn|,[£J. Moreover, the survival pro- 



3 If neutrinos are Majorana particles, there are n — 1 additional (so-called Majorana) phases, but they 
have no effect on neutrino oscillations. 
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babilities P aa , which are CP-symmetric in vacuum due to CPT invariance, are no longer 
CP-symmetric in matter. 

In contrast to this, the presence of matter does not necessarily break T invariance. In 
particular, matter of constant density (or, in general, matter with a density profile that 
is symmetric with respect to the interchange of the positions of the neutrino source and 
detector) does not induce any T violation. In addition, the survival probabilities P aa are 
always T-symmetric, since the initial and final neutrino flavours coincide. 

The point that symmetric matter density profiles do not induce any T violation seems 
intuitively rather obvious; here we will give an explicit proof of this statement and also will 
derive several useful properties of the evolution matrix S(t,t ). To study these properties, 
it is convenient to consider the evolution over the symmetric time interval [—t, t}. This does 
not lead to any loss of generality as any time interval can be reduced to a symmetric one by 
the proper choice of the point t = 0. Since both arguments of the evolution matrix S(t, —t) 
depend on t, it does not satisfy Eq. (|J). It is, however, not difficult to derive the evolution 
equation for S(t,—t). Using Eq. @ one can write S(t,—t) = S(t, 0)S(0, — t). The matrix 
S(t, 0) satisfies the usual evolution equation 

ij t S(t,Q) = H(t)S(t,Q). (9) 
Taking its Hermitian conjugate and substituting —t for t, one finds 

ils(0,-t)=S(0,-t)H(-t). (10) 
at 

Using Eqs. (^) and flHf), one finally obtains 

i-fS(t, -t) = H(t)S(t, -t) + S(t, -t)H(-t) . (11) 

The action of time reversal on this equation is given by the substitution t — > —t. 

Let us now assume that the fundamental CP and T violation is absent, i. e., all {Sep} = 0. 
In this case the Hamiltonian of the neutrino system is real (or can be made real by a 
rephasing of the neutrino states). Since it is real and Hermitian, it is also symmetric: 
H T = H. Transposition of Eq. (|ll|) yields 

ijS{t,-t) T = S(t,-t) T H(t) T + H(-tfS(t,-tf 

= S(t,-t) T H(t) + H(-t)S(t,-t) T . (12) 

The evolution matrix for neutrinos passing through the reverse profile N(t) = N(—t) 
can be obtained from Eq. (|ll~l) by replacing H(t) with H(—t). Comparing the resulting 
equation with Eq. (|T2"|), one obtains 

S&(t, t ) T = S iey (t, t ) ({5 CP } = 0) . (13) 
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Here the subscripts dir and rev denote propagation in the direct and reverse profiles, N(t) 
and N(t), respectively, and we have reinstated the time interval [t , t] in the arguments of the 
evolution matrices. From Eq. flHf ) it immediately follows that in the case of a symmetric 
matter density profile the evolution matrix S(t,to) is symmetric and there is no matter- 
induced T violation. For the particular case of matter consisting of a number of constant 
density layers, Eq. fll3|) was derived by Fishbane and Kaus [p3fl . 

It is easy to generalize Eq. (|13|) to the case when {Sep} 7^ 0. Following the lines that 
led to Eq. (O), one obtains, in this case, 



Sdir(Mo) — ^rev(Mo)|{ ( 5 Cp }.^_{ ( 5 Cp } • (14) 

This expression generalizes Eq. (13) of Ref. |23fl , which was obtained for matter consisting 
of constant-density layers, to the arbitrary matter density profile. Eq. flT4]) has a rather 
obvious physical meaning. It just reflects the fact that there are two kinds of effects that 
contribute to T violation (i.e., to the difference between S and S T ) in matter - intrinsic 
T violation, due to the non- vanishing CP and T violating phases {Sep}, and extrinsic T 
violation, due to the asymmetry of the density profile with respect to the interchange of the 
positions of the neutrino source and detector. 

Eq. (|13D means that in the case of {Sep} = the difference 

AP ab = P Aiv {v a -> u h ) - P vcv (u b -> v a ) (15) 

vanishes. Any deviation of this difference from zero is therefore a measure of non-vanishing 
fundamental CP and T violation and can, in principle, be used for their experimental 
searches. Fishbane and Kaus |23| have stressed that one can, in principle, probe the effects 



of {Sep} 7^ even by studying the survival probabilities P aa (a ^ e) if one compares these 
probabilities for direct and reverse profiles. This is a very interesting observation, even 
though the experiments with interchanged positions of neutrino source and detector would 
certainly be difficult to perform. 

The point that the survival probabilities can be used for looking for fundamental CP 
and T violation is easy to understand. The probabilities P{y a — > Va) are T-symmetric 
as the initial and final neutrino flavours coincide. Therefore, the contributions to their 
T asymmetry coming from the fundamental T violation and from the asymmetry of the 
matter density profile must cancel each other exactly. This means that by measuring the 
asymmetry P ( jir( z/ a —* v a) — Previa — ¥ v a) (with a ^ e) one directly measures, up to the 
sign, the asymmetry due to the fundamental CP and T-violating phases {Sep}- 

The asymmetry of the v e survival probability P ee with respect to density profile reversal 
vanishes (up to the tiny terms due to radiative corrections to matter-induced potentials of 

and v T [Q) due to the specific way the neutrino evolution equation depends on matter 
density p3| . Indeed, it was shown in Refs. fi, 25fl that P ee is independent of the phase 



Sep', from Eq. ( [HD it then immediately follows that P{y e — > ^ e )dir = P{y e — > u e ) Iev . It is 
interesting to note that this result depends crucially on the assumption that the matter- 
induced potentials of all neutrino species except u e are the same. This is no longer true if 
there are sterile neutrinos, and so in that case in general P{y e v e )&iv 7^ Pfae Z7 e)rev 
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An important point to notice is that there is no T violation (either fundamental or 
matter-induced) in two-flavour neutrino systems. Mathematically, this follows from the fact 
that the off-diagonal elements of any 2x2 unitary matrix have the same absolute values - 
for this reason the probabilities of the transitions v a —>■ uj, and z/f, — > v a , which are given by 
the squares of the moduli of the elements S21 and Si 2 of the evolution matrix S, coincide. 
Physically, this is related to the fact that conservation of probability (i.e., unitarity) puts 
rather stringent constraints in the two-flavour neutrino case. For example, in the (z/ e , z/ M ) 
system the conditions that the probabilities of the transitions from v e to all states (including 
v e itself) and from all states to v e must both be equal to unity are 

P ee + P eil = l, P ee + P tie = l. (16) 

From this one immediately obtains that 

i.e., neutrino oscillations are T- invariant irrespective of whether they take place in vacuum 
or in matter and whether the matter density profile is symmetric or not. 

Another consequence of unitarity in the two-flavour neutrino systems is that in the case 
of a symmetric matter density profile the off-diagonal elements of the evolution matrix S 
are pure imaginary. This is well known in the case of constant matter density, and it is easy 
to see that this in fact also holds for an arbitrary symmetric matter density profile. Indeed, 
in the two-flavour case the most general form of the unitary evolution matrix S (up to phase 
factors which can be absorbed into redefinitions of the fields) is 

S = ( £ ) , where |a| 2 + |/5| 2 = 1 . (18) 

From Eq. (|T3|) it follows that when N(t) = N(t), the evolution matrix is symmetric, i.e., f3 
is pure imaginary. 

Unitarity is much less constraining in the case of more than two neutrino flavours. For 
example, for three neutrino flavours one obtains from equalities similar to Eq. ( |TE| ) only the 
condition 

APj M + AP e T T = 0, (19) 

and the T-odd asymmetries AP e ^ and AP^ need not vanish. Considering different initial 
neutrino states, one can also find AP^ + AP^ = and AP^ + AP^ = 0, which together 



with Eq. (|T|) give § 

AP^ = APj T = AP T T e . (20) 

This relation coincides with the corresponding relation in the case of neutrino oscillations 
in vacuum; it implies that in a three-flavour neutrino system there is only one independent 
T-odd probability difference for neutrinos (and similarly one for antineutrinos). In a few 
recent publications |T5, 20"], this relation was obtained (in some approximations) for neutrino 
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oscillations in matter of constant density. However, as should be clear from its derivation, 
it is exact and does not depend on the matter density profile. 

It is well known that in the limits of vanishing mixing angle #13 (i.e., vanishing element 
U e s of the lepton mixing matrix U) or vanishing mass squared difference Am^i, three-flavour 
neutrino oscillations effectively reduce to the two-flavour ones | 26fl . Since, as was discussed 
above, there are no T-violating effects in the two-flavour neutrino case, any such effect can 
be considered as a measure of the genuine three-flavourness of the neutrino system. Studying 
T-violating effects in neutrino oscillations in matter can thus, in principle, provide us with 
important complementary means of measuring the parameters #13 and Am^, even in the 
absence of fundamental CP and T violation. 



3 Approximate analytic description for an arbitrary 
matter density profile 

We will give here a simple approximate expression for T-violating effects in neutrino oscilla- 
tions in a matter with an arbitrary matter density profile. This will enable us to study the 
interplay of the fundamental and matter-induced T-violating effects in neutrino oscillations. 

We will consider a three-flavour neutrino system [y e , v^, u T ) and use the parameterization 
of the leptonic mixing matrix U which coincides with the standard parameterization of the 
quark mixing matrix |[27||. The leptonic mixing angle #13 is constrained by the CHOOZ 



reactor neutrino experiment and is known to be small [28 



sin 2 2^3 <0.10 sin0i3<O.16. 



Analyses of solar and atmospheric neutrino data |29] also show that the ratio Am|/Am 



2 

atm 



AmyAmjj < 0.1. As was discussed above, T-violating effects (both fundamental and 
matter-induced) disappear in the limit of vanishing # 13 or Am|i and therefore they must be 
suppressed by both these small factors. This means that they can be calculated to a very 
good accuracy within perturbation theory. One can either consider Am^/Amjj as a small 
parameter and treat #13 exactly, or vice versa; the T-odd asymmetry will automatically 
be suppressed by both these parameters. We choose to treat sin# 13 as a small parameter, 
while making no further approximations. The details of our calculations are described in 
Appendix A; here we just describe the main idea and the results. 

In the standard parameterization, the leptonic mixing matrix U can be written in the 
form U = O23V13O12, where Oij are orthogonal matrices describing rotations by the angles 9ij 
in the corresponding -planes, and V13 is the unitary matrix which describes the rotation 
in the (l,3)-plane and in addition includes the Dirac-type CP and T-violating phase 5cp 
[see Eq. (|A3|)]. In the following, we will use the notation Sjj = sin^-, c^- = cos 9^, 

^^f 1 , A(t)EA-l[i + V(i)]. (21) 
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It is convenient to perform the rotation according to v' = Oj 3 z/, where v is the neutrino 
vector of state in flavour basis. To zeroth order in perturbation theory in the small parameter 
S13 = sin #13, the evolution matrix in the rotated basis can be written as 

/ a(t,t Q ) [3{t, t ) \ 
S (t,t )=[ -(3*(t,t ) a*(t,t ) , \a\ 2 + \P\ 2 = l, (22) 
V /(t,to) / 

where t 

/(*,*„) = expj-ijf A(f)df} (23) 

and the parameters a(t,t ) and /3(t,to) are to be determined from the solutions of the two- 
flavour neutrino problem in the (l,2)-subsector. It is now easy to obtain the correction to the 
evolution matrix S' to order s 13 , from which, upon rotation back to the unprimed basis, one 
finds the following expression for the T-odd probability difference AP e ^ = AP^ T = AP^ e []: 

API ~ -2 Sl3 s 23 c 23 (A - s\ 2 5) Im [e- iS °r/3*(A a - C* a )] , (24) 

where t t 

A a = a ! a*fdt' + (3 I (3*fdf, C a = f I af* dt' . (25) 

J to J to J to 

Eqs. (^4j) and ( |25|) give the T-odd probability difference in the three-flavour neutrino system 
in terms of / and the two-flavour neutrino amplitudes a and /3 for an arbitrary matter 
density profile. Note that this simplifies the problem considerably as the two-flavour neutrino 
problems are generally much easier to solve than the three-flavour ones. For 5c p = it can 
be shown that the right-hand side of Eq. (p4] ) vanishes for any symmetric matter density 
profile, as it must (see Appendix A). 

We will now consider two special cases: first, matter consisting of two layers of constant 
densities, and second, the adiabatic approximation for an arbitrary matter density profile. 
For the first case with two layers of constant densities we will use the following notation: 
layer widths L\ and L2, electron number densities Ni and A^2, the corresponding matter- 
induced potentials V\ and V2, the values of the mixing angle in the (l,2)-subsector in matter 
6\ and 9 2 0, respectively, 

Ui = ^(cos 29 12 5 - Vi) 2 + sin 2 29 12 5 2 , (26) 

and 

Si = sin(^Lj) , q = cos^UiLi) , (i = 1, 2). (27) 

4 Eqs. (Q) and ( ^5|) were obtained neglecting the (very small) corrections of the order of (Am^/Am^) 2 . 
Expressions which do not use this approximation are given in Appendix A. 
5 In the limit of vanishing densities N\ and N%, one has 6\ = 62 = Q\<x- 
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For the T-odd difference of the neutrino oscillation probabilities one then finds 
API ~ -2 Sl3S23 c 23 (A - s\ 2 5) [cos5 CP (X 1 Z R - X 2 Zj) + sin5 CP (X 1 Z 7 + X 2 Z R ) 



eft 



Here X\, X 2 , X 3 , and Y are the parameters which define the two-flavour evolution 
for neutrinos passing through the two-layer medium, s = Y — icr ■ X, with j30 



(28) 
matrix 



and 



Y = c\c 2 — S1S2 cos(20i — 26 2 ) , X\ = S1C2 sin 20i + S2C1 sin 29 2 , 
X 2 = -s 1 s 2 sin(26 l i - 28 2 ) , X 3 = -(sic 2 cos 29 1 + s 2 Ci cos 29 2 ) , 



= Re (A a - C* a ) = - [D_ - (fij cos 29 1 - Q 2 cos 2B 2 )\ 

x [sin(AiLi + A2-L2) - c 2 sin(AiL 1 ) - ci sin(A 2 -L 2 )] 

- [D_ cos 20 x + fi 2 cos(20i - 26 2 ) - ft^sj cos(A 2 L 2 ) 

- [£>__ cos 20 2 - fii cos(2#i - 20 2 ) + ^ 2 ]s 2 cos(AiLi) 

- L>_X 3 -fii[siC2 + s 2 cicos(20i - 20 2 )] + ^ 2 [cis 2 + sic 2 cos(20i - 20 2 )] 



, (29) 



Z 7 = Im (A a - C*) = [D+ - (fii cos 2^x + fi 2 cos 20 2 )] cos(AiLi + A 2 L 2 ) 

- [D_ - (Oi cos 26>i - fi 2 cos 29 2 )} [a cos(A 2 L 2 ) - c 2 cos(A 1 L 1 )] 

- [D_ cos 20i + Vl 2 cos(20i - 20 2 ) - Sl^si sin(A 2 L 2 ) 
+ [Z>_ cos 202-^1 cos (20x -20 2 ) + fi 2 ]s 2 sin(A 1 L 1 ) 

- D + Y + Qi(c\C2 cos 20i — sis 2 cos20 2 ) + fi 2 (cic 2 cos20 2 — sis 2 cos20i' 



(30) 



Here we have used the notation 



A, = A-i(i + V0, D ± = -^±^, a = ^ (i = l,2). (31) 



The term proportional to cos^cp in Eq. ( p8|) describes the matter-induced T violation, 
whereas the sin^cp term is due to the fundamental T violation. It is easy to see that the 
parameters X 2 and Zp are antisymmetric with respect to the interchange of the two layers, 
while Xx, X 3 , Y, and Zj are symmetric, and therefore the condition (|T4|) is satisfied. 

In the low energy regime 5 = Am 21 /2E > V 12 (which for the LMA-MSW solution of 
the solar neutrino problem and matter densities typical for the upper mantle of the earth 
corresponds to E < 1 GeV) the main contributions to Eq. (^) come from the D + terms in 
Zj, and the expression for the T-odd probability difference simplifies significantly: 

sinf20 20 ) 

APj ~ cos^cp ■ 8si 2 Ci2Si 3 s 23 C2 3 —{siS 2 [Y - cos(AiLi + A 2 L 2 )]} 

Sin /C7l2 

+ sin 5 CP • 4si 3 s 23 c 23 X 1 [Y - cos( A X L X + A 2 L 2 )] . (32) 
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Here the cos<5cp term has a remarkably simple structure: it is given by an oscillating term 
multiplied by an effective Jarlskog invariant 



_ sin(2gi - 29 2 ) 

sm 2&12 

This has to be compared with the usual Jarlskog invariant |yj 

J = Si2C 12 Si3Ci3S23C23 sin 5 C p . (34) 

Note that the factor sin(2$i — 29 2 )/ sin2#i2 in J c g plays the same role as sin Sep in J (which 
is zero in the absence of fundamental CP and T violation). The factor sin(2$i — 29 2 ) / sin 2# 12 
is a measure of the asymmetry of the matter density profile (6*i 7^ 9 2 ). Note also the absence 
of the c 2 3 factor which is equal to unity in our approximation. The sin 5c p term in Eq. ([32] ) 
can also be expressed in terms of the (usual) Jarlskog invariant if one writes X\ as 

y ■ 00 ( sin2g i sin2^2 \ 

X 1 = sm 26>i2 S1C2 . nn h s 2 ci . . (35) 

\ sm 2&12 sm 2012 / 

Note that the ratios sin2$i/ sin20i 2 and sin20 2 / sin20i 2 are finite in the limit sin20i 2 — * 0. 

Let us now consider an arbitrary matter density profile in the adiabatic approximation 
(see Appendix B). In this case the parameters a and (3 describing the two-flavour neutrino 
evolution in the (l,2)-subsector are 

a(t,t ) = cos$cos(6' - 6 ) + zsin$cos(0 + 9 ) , (36) 
P(t,t Q ) = cos$sin(0-0 o ) -zsin$sin(0 + o ), (37) 

whereas f(t,t ) is given by Eq. ( PB| ) as before. Here 8 and 9 are the values of the mixing 
angle in the (l,2)-subsector in matter at the initial and final points of the neutrino evolution, 
to and t, respectively, 

$= / u(t')dt', (38) 

•J to 

and uj(t) is given by Eq. (|26|) with the substitution Vi j2 — > V(t). We will also assume that 
V(t) < (5 <C A. The integrals in Eq. ( p5|) can then be done approximately (see Appendix B) 
and one obtains 

AP e ^ ~ 2S13S23C23 {cos 5 C p[sm(29 - 29 ) cos 2 $ - 2 sin(6» - 9 ) cos $ cos(A(t - t ))] + 

sin 5 C p[sm(9 + 9 ) cos(0 - 9 Q ) sin 2$ - 2 sin(0 + 9 ) sin $ cos(A(t - t ))]} . (39) 

In the regime in which the oscillations governed by large A = Amf 1 /2E are fast and 
therefore can be averaged over, the above expression simplifies to 

AP^ ~ 2S13S23C23 {cos 5 C p sin(20 - 29 ) cos 2 $ + sin 6 C p sin(0 + 9 ) cos(6» - 9 ) sin 2$} . (40) 
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Note that the cos^cp contribution can again be written in terms of the effective Jarlskog 
invariant J e ^\ 

sinf2$ 29 ) 

(APj) cos5cp ~ 4 cos 5 C p • S12C12S13S23C23 r^TT — — cos 2 $ = 4 cos 5 C p J eff cos 2 $ . (41) 

Sill ZU12 

Interestingly, the effective Jarlskog invariant appears both in the adiabatic regime and in the 
case of the two-layer matter density profile, which is an example of an extreme non-adiabatic 
case. If the oscillations due to a smaller energy difference u are also in the averaging regime, 
the sin 6cp contribution in Eq. fl40|) vanishes and only the matter-induced T- violating term 
proportional to cos^cp survives. Thus, we make the interesting observation that matter- 
induced T violation, unlike the fundamental one, does not disappear when the neutrino 
oscillations are in the regime of complete averaging. 



4 Implications 

In this section we compare our approximate analytic formulas with numerical calculations 
of AP e ^ and discuss the relevance of matter-induced T violation in neutrino oscillations for 
long baseline experiments, and also for solar, atmospheric and supernova neutrinos and for 
neutrinos in the early universe. 



4.1 Accuracy of the analytic approximation 

In order to estimate the accuracy and domain of validity of our results we have to distinguish 
two cases: L/E < 10,000 km/GeV and L/E > 10,000 km/GeV. In the first case the 
oscillating structure of APj b can, in general, be resolved, as it is shown in Fig. [!]. As can be 
seen from the figure, the size of AP^ is reproduced rather accurately, but there is an error 
in the oscillation phase which accumulates with distance and increases with increasing #13 
and Amli. 

In the second case, the oscillations governed by the large Am^ = Am 2 tm are very fast. 
We illustrate this case in Fig. |2|, where AP e ^ is plotted as a function of the total distance 
L traveled by neutrinos for matter consisting of two layers with densities p\ = 0, P2 = 6.4 
g/cm 3 and widths L\ = L 2 = L/2. We have chosen for the left plot the same parameters as 
those in Fig. 3a of Ref. p3| (]. The baseline values are rather unrealistic and serve illustrative 
purposes only. In the right plot of Fig. |2] we plot AP^ for larger values of # 13 and Am^. 



The result obtained using our analytic formula (p8|) (grey curve in Fig. |2|) reproduces the 



results of numerical calculations of Ref. 2^] and of our own numerical calculations very well: 
the difference can be barely seen. In fact, the fast oscillations cannot be resolved by any 



6 In Ref. |23| the value p 2 = 8 g/cm 3 instead of 6.4 g/cm 3 was erroneously quoted for Fig. 3a, but 
the quoted value of the potential V2 was correct. Note that in Fig. 3a of Ref. |E3[ the asymmetry of the 
transition probability P^ e with respect to density profile reversal rather than AP^ was plotted. However, 
as follows from Eq. (|13|), in the case of 5c p = this asymmetry coincides with AP^ e — —AP^ . 
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realistic detector due to its limited energy resolution; for this reason in Fig. |2| we compare 
the averaged over these fast oscillations values of AP^f obtained using approximate analytic 
formula (|28| ) with those calculated numerically (black solid and dashed curves, respectively). 
Since the error in the phase is irrelevant due to the averaging, the accuracy of the analytic 
approach is very good. For #13 = 0.1 and Am^ = 5 ■ 10~ 5 eV 2 the maximal relative error is 
about 5%; it increases with increasing # 13 and Am^, and for 13 = 0.16 and Am^ = 2- 10~ 4 
eV 2 the maximal error is about 10%. 

In both cases the accuracy of the predicted size of matter-induced T violation is quite 
good, the error being < 10%; however, the phase is not reproduced to the same accuracy. 




Figure 1: Comparison of the results of the analytic (solid curve) and numerical (dashed curve) 
calculations of AP^ for two layers of widths L\ = L2 = L/2, densities 1 g/cm 3 and 3 g/cm 3 and 
electron number fractions {Y e )\ = (Ye)2 = 0.5. In both plots E = 1 GeV, Am^ = 3.5 • 10 -3 
eV 2 , 5cp = 0, 6*12 = 0.56, and #23 = tt/4 are used. The remaining parameters are 6*13 = 0.1 and 
Amlx = 5 • 10 -5 eV 2 for the left plot and 6 13 = 0.16 and Am^ = 2 • 10~ 4 eV 2 for the right plot. 



4.2 Long baseline experiments 

One possible implication of the discussed effects may be in future long baseline neutrino os- 
cillation experiments. Neutrino factory experiments with thousands of kilometers baselines 
appear feasible, leading to beams which traverse the matter of the earth. The matter density 
profile of the earth is approximately spherically symmetric; however, the deviations from 
exact spherical symmetry may be considerable, which in turn can give rise to T-violating 
matter density profiles seen by neutrino beams. These fine details of the earth structure are, 
unfortunately, not very well known, but from seismological data one can infer some upper 
bounds on the possible asymmetry. A conservative estimate of these deviations is to assume 
a 10% variation in density on the length scales of several thousand kilometers |HJ some- 
where along the neutrino path. Another possibility to have an asymmetric matter density 
profile exists for shorter baselines L < 1, 000 km where the neutrino path goes only some 
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Figure 2: Comparison of the analytic and numerical calculations of AP^ for two layers of widths 
L\ = L2 = L/2, densities and 6.4 g/cm 3 and electron number fractions (Y e )i = (1^)2 = 0.5. The 
grey curve is the result of our analytic calculation, the black solid and dashed curves show the 
results averaged over the fast oscillations of the analytic and numerical calculation, respectively. 
In both plots E = 0.5 GeV, Am 2 31 = 3.5 • 10~ 3 eV 2 , 5 C p = 0, 6 12 = 0.56, and 6> 23 = vr/4 are used. 
The remaining parameters are #13 = 0.1 and Am|i = 5 • 10 -5 eV 2 for the left plot and 6*13 = 0.16 
and Amj! = 2 • 10 -4 eV 2 for the right plot. 



ten kilometers below the surface of the earth. In this case one half of the neutrino path 
could be in the sea (p — 1 g/cm 3 ) and the other half in the continental earth crust (p — 3 
g/cm 3 ). The experimental setups under discussion will have a remarkable precision, and it 
is thus important to know when the discussed matter-induced T-violating effects become 
relevant. In other words, we will discuss when it is necessary to include the discussed extra 
T-violating effects and when the usual analysis is justified. 

We consider the usual setup for neutrino factory experiments with a beam energy of 50 
GeV. The energy threshold of the detector is 4 GeV and its energy resolution is 10%. The 
detector is capable of charge identification, thus the wrong sign muon signal is available. 
We assume 2 ■ 10 21 useful muon decays for both polarities and a detector mass of 40 kton. 
This luminosity coincides with the one used in Ref. |I| and is a factor of 40 larger than the 
one used in Ref. 0. We include only statistical errors. All fits include both polarities and 
both appearance and disappearance rates. Further details can be found in Ref. 0. 

We assume two different asymmetric matter density profiles. The first one consists of 
two layers of equal widths, with densities 1 g/cm 3 and 3 g/cm 3 (Figs. |3| and f|). This 
corresponds to the sea-earth scenario, which of course can be realized on the earth only up 
to baselines of ~ 1, 000 km. The second one also consists of two layers of equal widths, but 
with densities 3 g/cm 3 and 3.3 g/cm 3 as an example of density perturbations which could 
arise in real very long baseline experiments (Fig. ^|). 

We simulate such experiments numerically for these types of profiles and perform fits to 
the obtained event rates. We compare this with the fits performed for symmetrized versions 
of the corresponding profiles, which are modeled by replacing the transition probabilities by 



14 



the symmetrized ones: 

P S =\ (Pdir + Pre.) ■ (42) 

Thus we are only sensitive to the errors induced by the asymmetry of the matter density 
profile and not to the errors in the average density, which of course should not be ignored in 
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Figure 3: A fit of sin 2 2^13 and 8c p for different baselines in a neutrino factory experiment with 
a beam energy of 50 GeV (see the text for details). The solid curve is the 99% C.L. contour for 
the asymmetric matter density profile with two layers of widths L\ = L2 = L/2, densities 1 g/cm 3 
and 3 g/cm 3 and electron number fractions (Y e )i = (1^)2 = 0.5. The dashed curve is the 99% C.L. 
contour for the symmetrized matter density profile as defined in Eq. (|42|). The star and square 
indicate the best fit points for the asymmetric and symmetrized profiles, respectively. The Ax 2 
value given in each plot is the difference of x 2 (2 d.o.f.) between the two best fit points. The 

12 = 7T/4, 6*23 

and 5cp = 0. 



parameters are d 12 = tt/4, 6*23 = vr/4, Am§ 2 = 3.5 • 10" 3 eV 2 , Aro^ = 10 -4 eV 2 , sin 2 29 13 = 0.1 
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Figure 4: Same as in Fig. |||, but for the fixed baseline L = 6, 000 km and three different values of 
5cp- — 7r/2, 0, and ir/2. 
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reality. The difference between the minimal \ 2 values for the asymmetric and symmetrized 
profiles is a direct measure of the sensitivity to the matter-induced T violation. 




The effect can be quite sizable for matter density profiles of the sea-earth type (see 
Fig. |3p, however, only for baselines above 1,000 km, which cannot be realized on the earth. 
For large baselines the errors in the determination of the fundamental CP and T-violating 
phase 5c p induced by asymmetric matter are comparable with the statistical errors in the 
case of symmetric matter. This behavior is quite similar for all possible values of 5c p as 
can be seen in Fig. |j. 
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Figure 6: Energy dependence of the transition probability P efM (grey dashed curve), T-odd proba- 
bility difference AP e ^ multiplied by 10 (grey solid curve), and the ratio AP^/(P efM ) , where (P ef i) 
is the average of P e ^ over fast oscillations (black solid curve) in the case of two-layer density profile 
with pi = 1 g/cm 3 , p2 = 3 g/cm 3 and electron number fractions (Y e )i = (Y e )2 = 0.5. The layer 
widths are L\ = L<i = 500 km (left plot) and L\ = L2 = 1,500 km (right plot). The values of 
neutrino parameters are Am^ = 2 • 10~ 4 eV 2 , #12 = 0.559, Am^ = 3.5 • 10~ 3 eV 2 , #23 = vr/4, 
6»i3 = 0.16, and 5 C p = 0. 
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The energy dependence of the transition probability P e(U , T-odd probability difference 
AP e ^, and their ratio in the case of two-layer matter density profile corresponding to the 
sea-earth scenario is shown in Fig. |[ One can see that the relative size of the T-violating 
effect is largest at energies of about 1 GeV. A dedicated low energy experiment capable of 
measuring the matter-induced T violation with the sea-earth type matter density profile at 
a baseline of ~ 1, 000 km would, however, require enormous luminosities, at least six orders 
of magnitudes higher than an initial stage of a neutrino factory (2 ■ 10 19 muons/year, 10 
kton) . 

For much more realistic matter density profiles with only 10% density variation, the 
matter- induced T violation effects are small at any baseline (see Fig. [5]). The statistical 
errors in the determination of 8c p and #13 are much larger than the errors induced by 
replacing the asymmetric matter density profile by the symmetrized one. This again holds 
for all values of 8c p- Therefore we conclude that the determination of the fundamental CP 
and T-violating phase 8c p cannot be spoiled by an unknown T asymmetry of the matter 
density profile until the luminosity (and therewith the accuracy) is increased by at least two 
orders of magnitude, bringing the error in 8c p down to ~ 1% at 99% C.L. 

4.3 Solar, supernova and atmospheric neutrinos and neutrinos in 
the early universe 

Matter-induced T-violating effects can, in principle, manifest themselves in a number of 
other situations when oscillating neutrinos propagate through asymmetric matter. In par- 
ticular, solar neutrinos traveling from the center of the sun towards its surface traverse a 
matter density profile which is highly asymmetric and characterized by a large contrast of 
densities (from p ~ 150 g/cm 3 to p ~ 0). There are no muon or tau neutrinos origi- 
nally produced in the nuclear reactions in the sun, and so one cannot study the quantities 
like AP^ or AP^ T with solar neutrinos f\ However, genuine three-flavour effects (including 
those of asymmetric matter) in general contribute terms of the order of S13 to the transition 
probabilities which can be large compared to the usual terms proportional to s^ 3 which are 
present in the effective two-flavour approach. In particular, the survival probability of solar 
neutrinos averaged over the fast oscillations due to the large Am^ tm is approximately given 
by "quasi two-flavour" formula |BB 



Pee = 4 3 P£ + St 3 , (43) 

where is the two-flavour v e survival probability in the (1,2) sector calculated with the 
effective matter- induced neutrino potential V e s = cl 3 V. It contains only second and higher 
order in S13 corrections to the standard two-flavour result. One can then ask if this result 
is modified by asymmetric matter, leading to the order S13 corrections. 

7 Neutrinos of all flavours with energies O(GeV) can be produced in the annihilation of weakly interacting 
massive particles (WIMPs) inside the sun. T-violating effects of solar matter on the oscillations of such 
neutrinos were discussed in Ref . I[lq| . 
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In the case of the LMA-MSW solution of the solar neutrino problem (which leads to the 
largest allowed values of Am^) the evolution of neutrinos inside the sun is adiabatic, and 
Eqs. (|3^) or fl40|) can be used. For maximal values of #13 allowed by the CHOOZ experiment 
one finds that the corrections of the order S13 to the transition probabilities P efl and P eT 
of solar neutrinos can be as large as (5 - 10)%. However, these corrections have opposite 
sign and the same absolute value, and therefore they cancel exactly in the total transition 
(or survival) probability of electron neutrinos. The fact that the survival probability of 
v e is symmetric with respect to matter density profile reversal has already been discussed 
above (see Sec. 0); since one cannot distinguish from v T in low energy experiments, 
this makes the asymmetric matter effects in oscillations of solar neutrinos inside the sun 
unobservable. Eq. (pE3|) therefore does not have to be modified. We have checked this by 
comparing its predictions with the results of numerical calculations of the full three-flavour 
evolution equation for solar neutrinos and found that the accuracy of the approximation in 
Eq. ([43|) is extremely good (error < 10~ 3 ). 

Another potentially interesting implication of matter-induced T-violating effects could 
be in oscillations of supernova neutrinos propagating from the supernova's core outwards. 
However, these effects again cancel because of the absence of the matter-induced asymmetry 
of P ee and the fact that the fluxes and spectra of the supernova and v T are practically 
identical. 

The expansion of the universe implies that neutrino oscillations in the early universe take 
place in a time-dependent, asymmetric environment. One can then ask if T violation effects 
in neutrino oscillations can result in a lepton flavour asymmetry. Such an asymmetry, e.g., 
in electron neutrinos could have important consequences for the big bang nucleosynthesis. 
However, it is easy to make sure that if the original numbers of neutrinos of all species are 
equal, this does not happen: the summation over the contributions of all flavours makes the 
asymmetries in the individual lepton flavours vanish. 

Thus, we see that even though there are, in principle, interesting extra T-violating effects 
in asymmetric situations like in the sun, in the early universe and in supernovae explosions, 
these effects cancel whenever a summation over neutrino flavours is present, and a naive 
analysis without these effects leads to correct results. 

Matter-induced T violation can, in principle, also influence atmospheric neutrino oscil- 
lations. This is related to the fact that neutrinos are produced in the atmosphere at an 
average height of about 15 km and so the neutrinos coming to the detector from the lower 
hemisphere travel first in air and then in the earth. For neutrinos from just below the hori- 
zon this can result in nearly equal pathlengths (of the order of a few hundred kilometers) 
in air and in the earth. The effect is, however, very small, which can be seen through a 
simple estimate. The asymmetric matter effects can be looked for through the distortions 
of the zenith angle distributions of the e-like events (the distributions of the /i-like events 
are mainly governed by the large 623 and Am^, and the relative T-violating effects for them 
are smaller than those for the e-like events). The probability difference AP e ^ is of the order 
of a few percent of P eAt , and P efM itself is of the order of a few percent of P ee ; therefore the 
matter-induced T-violating effect is of the order 10~ 4 — 10~ 3 at best. Thus, one would need 
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10 6 — 10 8 events to achieve a statistical accuracy of the same order. Super-Kamiokande at 
the moment has roughly 10 4 events. Including the flux uncertainties would make the situa- 
tion even worse. Therefore matter-induced T- violating effects in atmospheric neutrinos can 
be safely neglected. This result is also confirmed by numerical analysis. 

5 Discussion and conclusions 

In the present paper we studied T violation in neutrino oscillations in asymmetric matter, 
with a special emphasis on matter-induced T violation. To this end, we derived a simple 
approximate analytic expression for T-odd probability differences AP^ in the case of a 
general matter density profile. We have shown that our analytic expressions reproduce the 
results of direct numerical integration of the neutrino evolution equation very well. 

Matter-induced T violation has two aspects. First, it is an interesting matter effect, 
which is present only in asymmetric matter (and therefore is absent in the most studied 
case of constant-density matter). It can manifest itself not only in the T-odd differences 
of oscillation probabilities but also in specific modification of the probabilities themselves. 
Since this effect can only exist in systems of three or more neutrino flavours, it is sensitive to 
the parameters that discriminate between genuine three-flavour and effectively two-flavour 
oscillations, such as #13 and Am^, and therefore can, in principle, be used for their deter- 
mination. Second, matter-induced T violation can fake the fundamental one and so impede 
the determination of the fundamental CP and T- violating phase 5c p in the long baseline 
experiments. 

We have studied both these aspects and have shown that in the case of three neutrino 
species and for matter density contrasts and baselines feasible in terrestrial experiments, 
matter-induced T violation effects are small and can safely be ignored. This is due to the fact 
that they are doubly suppressed by small factors sin #13 and Am^/Am^ = Am|,/Am^ tm . 
In particular, asymmetric matter effects cannot hinder the determination of the fundamental 
CP and T- violating phase 5c p in the long baseline experiments as far as the error in the 
determination of 5cp is larger than 1% at 99% C.L. This is the main result of our paper. 

T-asymmetric matter effects for solar, supernova, and cosmological neutrinos, which 
propagate through larger density contrasts and over larger distances, can, in principle, be 
large in each individual oscillation channel. However, these effects are not observable because 
of the summation over channels inherent in the experimental detection of these neutrinos 
or in the observable quantities such as total number density of neutrinos of a given flavour 
in the case of neutrinos in the early universe. 

The situation can be very different in four-neutrino schemes. In that case T-violating 
effects (both fundamental and matter-induced) are not, in general, suppressed by small 
factors like sin 13 or Am|/Am^ tm . Therefore experimental detection of sizeable T violation 
effects would signify the existence of a fourth light neutrino, which, due to the LEP result 
on the invisible Z° boson width, must be a sterile neutrino. It should also be understood 
that our conclusion that matter-induced T violation cannot hamper the determination of 
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the fundamental CP and T- violating phases does not, in general, apply to the four-neutrino 
schemes; in this case one has to carefully take matter-induced T-violation into account in 
order to disentangle intrinsic T violation from the extrinsic one. 

It is also interesting to note that the statement that the survival probability of electron 
neutrinos P ee is invariant with respect to matter density profile reversal, discussed in Sec. 0, 
holds only in three-neutrino schemes. Indeed, it relies on the assumption that the matter- 
induced potentials of all neutrino species except u e are the same. This is no longer true if 
sterile neutrinos are present. Therefore an observation of non-invariance of P ee under density 
profile reversal would be a signature of a sterile neutrino. It should be noted, however, that 
both this non-invariance of P ee and large T-violating effects can only be observable in the 
(experimentally favored) 2+2 four-neutrino schemes, in which there are two pairs of nearly 
degenerate neutrino mass eigenstates separated by a large mass gap. In the 3+1 scheme 
with a lone neutrino mass eigenstate being predominantly a sterile neutrino these effects 
are expected to be small since this scheme can be considered as a small perturbation of the 
three-neutrino ones. 
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Appendix A: Analytic description of AP^ 

The evolution equation describing neutrino oscillations in the three-flavour neutrino system 
can be written as 
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where u a (a = e, fi, r) are the components of the neutrino vector of state in flavour basis, U 
is the leptonic mixing matrix, 5 and A were defined in Eq. (pTJ), and 



V(t) = V2G F N e (t) 



(A2) 



is the charged-current contribution to the matter-induced potential of electron neutrinos, 
Gf and N e being the Fermi weak coupling constant and the electron number density of 
the medium, respectively. The neutral-current contributions to the potentials of v e , v. 



and v T in matter are the same (up to tiny radiative corrections 0]) and therefore do not 
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affect the neutrino oscillations probabilities. The evolution matrix S(t, to) satisfies the same 
Schrodinger equation ( |A1| ). 

In the standard parameterization, the leptonic mixing matrix U can be written in the 
forniQ U = O23V13O12, where 

/ ci2 % 0\ / cia s 13 e~ i5 ^ p \ / 1 

12 = s 12 cia , V 13 = I , 23 = c 23 s 23 

V 1/ V -S 13 e^ C13 / V " s 23 c 23 

(A3) 

Since the matrix of matter- induced neutrino potentials diag (V(t), 0, 0) commutes with 23 , 
it is convenient to perform the rotation according to v' = 0\ ? v (or S' = 23 S0 23 ), where 
v = {y e v T ) T and S are the neutrino vector of state and evolution matrix in flavour basis. 
The effective Hamiltonian H'(t) governing the neutrino evolution in the rotated basis can 
be written as 

/ cj 3 sj 2 5 + s 2 13 A + V(t) c 13 c 12 s 12 5 c 13 s 13 (A - s\ 2 5) e~ iS ^ \ 
H\t) = c 13 c 12 s 12 5 c 2 12 5 s 13 c 12 s 12 e- iS ^5 . (A4) 

\ C13S13 (A - s\ 2 5) e i5cp -s 13 c l2 s l2 e i5cp 5 s 2 n s\ 2 8 + cf 3 A / 

It can be decomposed as 

J ff'(t)=^(t) + ^, H'j = H' 1 + H^ (A5) 

where 

/ |[- cos 20ia <J + y(t)] 2 -sm29 12 6 \ / ( ) 

H' (t)=l ' \sm26 l2 5 |[cos 29 l2 5 - V(t)] = [ ) 

\ A(t) / \ A(t) 

(A6) 

is of zeroth order in the small parameter S13, H[ is of the first order, and H 2 includes terms 
of the second and higher orders, namely 

/ sis (A - s\ 2 8) e~ iScp \ / a 

H[=\ -si 3 ci 2 si 2 e-^5 =006 

\ si 3 (A - s\ 2 5) e l5cp s l3 c 12 s l2 e i5cp 5 / \ a* b* 

(AT) 

and H 2 = 0(s\ 3 ). The function A(t) that enters into Eq. ( |A6[ ) was defined in Eq. (PTT) . 
We have subtracted from the Hamiltonian a term proportional to the unit matrix in order 
to make the upper-left 2x2 submatrix h of H' Q traceless. This amounts to multiplying all 
the components of the neutrino state by the same phase factor, which does not affect the 
neutrino oscillation probabilities. 

We will look now for the solution to the Schrodinger equation for the evolution matrix 
5" in the rotated basis in the form 

S'(t,t ) = S' (t,t )S[(t,t ), (A8) 



^Omitting possible Majorana phases which have no effect on neutrino oscillations. 
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where S' (t,t ) satisfies the equation 

ij t S' (t,t ) = H' (t)S' (t,t ) , (A9) 

with the initial condition S' Q (to,to) = 1. Its general solution can be written as 

a(t,t ) P(t,to) 
S' (t,t )=\ -P*(t,t ) a*(t,t ) |, \a\ 2 + \/3\ 2 = l, (A10) 
f(t,to) 

where f(t, t ) is given in Eq. fl23|) and the parameters a(t, t ) and j3(t, t Q ) are to be found from 
the solution of the two-flavour neutrino problem that corresponds to the 2x2 submatrix h 
of H' . As follows from the evolution equation for S' and Eqs. ( |A8| j and ( |A9| ) , the matrix S[ 
satisfies the equation 

ij t S[(t,t ) = [Sfah^HiS&MS&M , S[(t ,t ) = 1. (All) 

Up to now everything is exact; we will determine now the evolution matrix S' to the first 
order in perturbation theory in the small parameter S13. This leads to 

S'(t,to) ^ S' (t,t )-iS' (t,t ) f[S' Q {t',Uy l H[S' {t',t Q )}dt' . (A12) 



to 



A straightforward calculation then gives 



a(Mo) /3(Mo) -iA 
S'(t,t )=\ -/3*(t,t ) a*(t,t ) -iB |, (A13) 
-iC -iD f(t,t ) 



where 



A = aA a + bA b , B = aB a + bB b , (A14) 
C = a*C a + b*C b , D = a*D a + b*D b , (A15) 

the parameters a and b were defined in Eq. (|A7|) and 

A a = a(t,t ) j a(t',t )*f(t',t )dt' + P(t,t ) j P{t',t Q )*f{t',t Q )dt', (A16) 

A b = p(t,t ) f a(t',t Q )f(t',t )dt' -a(Mo) / P(t',t )f(t',t )dt' , (A17) 

5 a = a(Mo)* / P{t',t yf{t',t )dt' -P{t,t )* [ a(t',t yf{t',t )dt' , (A18) 
-'to -'to 

5 6 = a{t,t )* I a{t',t )f{t',t )dt' + P{t,t o y I P{t',t )f{t',t )dt', (A19) 

./to -/to 
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as well as 



C a = f(t,t ) [ a(t',t )f(t',t )*dt', C b =-f(t,t )[ P(t',t yf(t',t )* dt', (A20) 

J t J t 

D a = /(Mo) / P(t',to)f(t',toydt', D b = f(t,t ) I a(t',t yf(t',t o ydt'. (A21) 



'to •'to 

The calculation of the integrals in Eqs. ( A16| )-( [Al9| ) can be simplified considerably by notic- 
ing that their right-hand sides contain the expressions that are the products of the elements 
of the evolution matrix S' (ti,t). Indeed, one has 

( a x Pi \ / a* -(3 \ 
S' (t u t) = S' (t u t )S' (t,t o y = \ -Pt at )[ 0* a , (A22) 

V o o h ) \ o o/*/ 

where we used the notation a = a(t,to), a± = a(ti,t ), etc. On the other hand, we have 

/ a(t u t) p(t u t) \ 
S' {t 1 ,t)=[ -P(h,ty a(h,ty . (A23) 
V o f(h,t) J 

Comparing Eqs. ( |A22| ) and (|A23|) one finds 

(aat + PPDfJ* = a{h,tyf(t 1} t) , (Pa, - ap 1 )f 1 f* = -P(h,t)f(t 1} t) , (A24) 
(a* PI - P*a{)f 1 f* = P(t u tyf(t u t) , (a*ai + P*Pt)f 1 f* = <*(*i, t)/(*i, t) . (A25) 

This allows a simplification of the integrals in Eqs. (|A16[) -( [Al9|) , making them similar in 
form to those in Eqs. ( |A2D|) and QA2T] ): 

A a = fit, t ) f a(t', t)*f(t', t) dt' , A b = -fit, t ) f P(t', t)f(t', t) dt' , (A26) 

J to Jt 
B a = -f(t,t ) j P(t',tyfit',t)dt', B b = fit,t ) I a(t',t)f(t',t)dt' . (A27) 

Jto Jto 

Rotating S'(t, i ) by 23 back to the original flavour basis, one finds from Eq. (|A13| ) 

(a c 23 P - IS23A -s 23 P - ic-xA \ 

-c 23 /3*-is 23 C S 22 S 23 , (A28) 

s 23/3* — 2C 23 C S32 S33 J 

where 

5 22 EE 4 J a* + 4,/-is 2 3C 23 (5 + £>), (A29 ) 

5 23 EE -S23C23(a*-/)-i(45-*23^)> ( A 30) 

5 32 EE -5^3(0;* _/)+i( S 2 3jB _ c 2 3jD)) (A31) 

5 33 ee 4a* + 4/ + is 23 c 23 (B + D). (A32) 
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The parameters a, ft, and / in Eqs. ( |A13| ) and (|A2S| )-( |A~32|) are of zeroth order in s 



'1.3 



whereas A, B, C, and D are of the first order. From Eq. ( |A13| ) one can see that, if the 
mixing angle 6*23 were zero, the T-odd probability difference AP^ e = \A\ 2 — \C\ 2 would have 
scaled with #13 as s 2 3 f\ We have checked this by solving Eq. (|A1|) numerically. Since in 
vacuum there is no CP and T violation when any of the mixing angles is equal to zero, 
non-vanishing APj b = \A\ 2 — \C\ 2 in the case #23 = is a pure matter effect. Moreover, it is 
an effect of asymmetric matter: as we show below, in the case of matter with a symmetric 
density profile (and so also in vacuum) \A\ 2 = \C\ 2 , and all APj b vanish. The atmospheric 
neutrino data indicate that #23 is close to 45°, and therefore the T-odd probability differences 
must scale linearly with S13. From Eq. (|A28 ) one finds 



AP e ; = l^xl 2 - \S 12 \ 2 = -2s 23 c 23 lm [ft*(A - C*)} . (A33) 

From the definition of the parameters a and b in Eq. ( |A7| ) it follows that the ratio \b/a\ ~ 
AmyAffijj is small; it can also be shown that | v4& / A a | ~ \Cb/C a \ ~ Am^/ Am 3l . There- 
fore, the contributions of A b and C b to Eq. ( |A33| ) are suppressed by the factor (Am^/Am^) 2 , 



and one finally arrives at Eq. (|24|). 

We will now show that our expression ( [A3 3]) satisfies the requirement that in the absence 
of the fundamental CP and T violation (i.e., in the case Sep = 0), the T-odd probability 
differences must vanish for any symmetric matter density profile. To do so we will show 
that in this case the expression ft* (A — C*) is real. 

Consider the evolution of the neutrino system over the symmetric time interval [—t, t] . 
First, from the evolution equations (|9J) and fllOD we notice that in the case of symmetric 
matter the matrix S' (0, — t) T coincides with S' (t, 0) (the phase 5cp does not affect the 
evolution equations for S' ). Together with Eq. (|8|) this leads to the following symmetry 
properties of a, ft, and / in symmetric matter: 

a(-t, 0) = ait, 0)* , ft(-t, 0) = ftit, 0)* , /(-*, 0) = f(t, 0)* . (A34) 

Consider now the evolution matrix S'(t,—t) = S'(t, 0)S"(0, — t). It can be parameterized 
in the form similar to that of Eq. ( |A13| ), and its entries can be found in terms of those of 
S'{t, 0) and S"(0, —t). Using Eq. ( |A34j) it is then straightforward to show that in symmetric 
matter the entries of S(t, —t) satisfy 

A a = C a , A b = C b . (A35) 



In the case of 5c p = the parameters a and b are real, and from Eqs. ( |A15| ) and ( |A35|) one 



finds A = C . This means that A — C* is pure imaginary; since ft is also pure imaginary in 
this case, the right-hand side of Eq. (|A33| ) vanishes, which completes the proof. 



As follows from the definition of the parameters a and b in Eq. ( |A7[ ), they can be written 
as a = e~ l5cp a' , b = e~ lScp b' with real a' and b'. Therefore in the general case of 5cp 7^ 
one finds that in symmetric matter \A\ = |C|. 



9 Duc to Eq. ( pot ) so would AP^ and AP^. do. However, the former cannot be found directly from 
Eq. ( Al3| ) as this would require the calculation to be done in the next order in perturbation theory. 
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Appendix B: Two special cases 



We shall now calculate AP e ^ in two special cases - matter consisting of two layers of constant 
densities and the adiabatic approximation in the case of an arbitrary matter density profile. 

Consider first matter consisting of two layers of constant electron number densities Ni 
and N2 and widths L\ and L2, respectively. The corresponding matter-induced potentials 
(|A2|) are V\ and V2. The values Q\ and 82 of the mixing angle in the (l,2)-subsector in matter 
of densities N% and N 2 are given by 

cos 26 12 5 - Vi cos29 12 5 -V 2 , . 

cos26'i = , cos 28 2 = , (Bl) 



where w\ and UO2 were defined in Eq. (p6|). The time interval of the evolution of the neutrino 
system can be divided into two parts: (I) < t' < L\ and (II) L\ < t' < L\ + L 2 = L. 
The parameters a, (3, and / for the first interval are given by the well-known evolution in 
matter of constant density: 

a(t', 0) = cos^it') + i cos 26 1 sin^t') , (3(1/, 0) = -i sin 26 1 sin(cj 1 t') , (B2) 

f(t', 0) = exp(-iAif) , Ai, 2 = A - ^(5 + V 1>2 ) , (B3) 

whereas for the second interval they can be expressed through the elements of the two-flavour 
neutrino evolution matrix in the two-layer medium [[30]: 

a(t', 0) = cic' 2 - sis 2 cos(26'i - 20 2 ) + i(sic' 2 cos 26 l + s 2 ci cos 29 2 ) , (B4) 

(3(t\ 0) = -i(sic 2 sin 29 x + s 2 ci sin 29 2 ) + s x sf 2 sin(20i - 20 2 ) , (B5) 

/(*', 0) = expl-zfAiLi + A 2 (t' - Li)]} . (B6) 

Here 

s 2 = sin(a;2r) , c 2 = cos(u;2t) , r = t' — L± , (B7) 

and 51,2 and ci^ were defined in Eq. ( p^) . Direct calculation using Eqs. (|2^), ( |A26| ) , and 
Q2|) [or Eqs. (g), (^T§), and (g§] gives the result presented in Eqs. (|2|)-(|3p]). 

Next, we consider the adiabatic regime in the case of an arbitrary matter density profile. 
In this regime, the change of the matter density along the neutrino path is slow compared to 
the oscillation frequency. The effective Hamiltonian of the neutrino system is approximately 
diagonal in the basis of instantaneous matter eigenstates, and the evolution of the eigenstates 
amounts to a mere multiplication by phase factors. The evolution in the flavour basis is 
obtained by rotating the evolution matrix in the instantaneous eigenstates basis S(t, to)ei g en 
by the matrices U(to) and U(t) of leptonic mixing in matter, which correspond to the 
initial and final times of neutrino evolution: S(t,t ) = U(t)S(t, to)eigen^ r (^o)^- Applying 
this procedure to the two-flavour neutrino evolution described by the Hamiltonian Hq, one 
arrives at Eqs. (|36| ) and (j3~?D , while the parameter / is always given by Eq. (p3|) . 

We will be assuming that V(t) < 5 A, and also that the oscillations governed by 
the large A = Am'^ l /2E are in the averaging regime (i.e., that A • L ^> 1). The latter 
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assumption allows us to obtain simple approximate expressions for the relevant integrals. 
We first note that 

jf g(t')e' iAt ' dt> ~ ^[g(t)e- lAt - g(t )e- iAt °) + 0(1/ A 2 ) , (B8) 

where g(x) is an arbitrary regular function which changes slowly on the time intervals ~ 1/A, 
and the integration by parts has been used. Direct calculation then yields 

A a - C* ~ z-|{cos[A(t - t )] - cos{6 - ) cos $} . (B9) 

Eqs. d|), (|36|), and QB?D lead to Eq. (||). 
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